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The Hidden Local Symmetry Lagrangian is used to study the interactions of (j) mesons with 
other pseudoscalar and vector mesons in a hadronic gas at finite temperature. We have found a 
significantly small (p mean free path (less than 2.4 fm at T > 170 MeV) due to large collision rates 
with p mesons, kaons and predominantly K* in spite of their heavy mass. This implies that </!> 
mesons produced after hadronization in relativistic heavy ion collisions will not leave the hadronic 
system without scattering. The effect of these interactions on the time evolution of the (p density in 
the expanding hadronic fireball is investigated. 
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PACS numbers: 25.75.-q, 14.40.Ev, 24.10.Pa, 12.39.Fe 



I. INTRODUCTION 

Enhanced strangeness production in relativistic 
nucleus-nucleus collisions compared to nucleon-nucleon 
collisions has been predicted to be a consequence of the 
formation of quark-gluon plasma (QGP) (l). This effect 
leads to an enhancement in the number of strange and 
multistrange particles produced after hadronization. In 
particular, free ss pairs would coalesce to form (j) mesons, 
while their production in pp and np collisions should be 
OZI suppressed 1^. On the other side, the (j) mass is 
expected to decrease due to both many body effects in 
the hadronic medium Q and chiral symmetry restora- 
tion, which might generate a double peak structure 
The decay width is also modified (broadened) through 
the scattering with other particles [|[ 0, ||. 

The is a nice probe since it is not masked behind 
other resonances in the mass spectra. It decays into kaon 
pairs {K~^K~), and more rarely into dileptons (e~'"e~, 
both channels have been detected at CERN- 
SPS H, |lO|, Dileptons have negligible final state in- 
teractions with the hadronic environment, so they sense 
the entire evolution of the system. On the contrary, de- 
tectable kaons from (j) decay probably emerge only at 
freezeout. 

It is widely accepted that the 4> mean free path (MFP) 
in the hot hadronic fireball is large due to the small 
cross section for scattering with nonstrange hadrons. 
This implies that (p spectra would retain the informa- 
tion about the stage of the collision at which the plasma 
hadronizes Available calculations D seem to 

support this idea. For example, in Ref. the (j) MFP 
was calculated taking into account its scattering with 
different mesons. Phenomenological Lagrangians with 
couplings extracted from the experimental partial decay 
rates were used. The obtained MFP at T 200 MeV is 
rather big (A = 4.4 fm) compared to the standard size of 
the hadronic system. Adding the reactions with nucleons 
and nucleonic resonances did not change qualitatively the 
situation ||]. 

However, (p production in Pb — Pb collisions at SPS 



shows some intriguing features that are difhcult to match 
with the picture of a </) weakly interacting with the 
hadronic medium. Both absolute yields and inverse slope 
parameters in the transverse mass (m^) distributions ex- 
hibit different values when measured via or K~ 
decays [|T^. The inverse slope, as obtained from the 
hadronic measurement, suggests that the 0's flow to- 
gether with pions, kaons and protons, while the dilep- 
ton measurement is consistent with the assumption of an 
early freezeout. It has also been observed that the rapid- 
ity distribution (extracted from kaon pairs) in Pb — Pb 
is about 50 % broader than in ppjst]. The origin could 
be attributed to longitudinal flow ||l^ . The modification 
of visible spectra due to kaon rescattering inside the 
fireball is an important ingredient but does not fully ex- 
plain the discrepancies. Indeed, Johnson et al. |0 show 
that kaon rescattering can account for the observed ra- 
pidity distributions, but not for the differences in the rrit 
spectra, neither the inverse slopes, nor the relative yields. 

Here we calculate (j) collision rates and mean free path 
in a hot hadronic gas of pseudoscalar (tt, K) and vector 
mesons (p, lo, K* , (f>). The reaction cross sections are 
obtained within the Hidden Local Symmetry Lagrangian 
(HSL) |l^, which includes both Goldstone bosons and 
vector mesons in a manner consistent with the symme- 
tries of QCD. The use of such a realistic model allows 
us to take into account many mechanisms that are not 
present in calculations that rely only in couplings ex- 
tracted from observed decays but are allowed by the sym- 
metries. This is, for instance, the case of vertices like 
(j>K*K, pK*K* and many others. As a consequence, we 
shall see that at temperatures between 150 and 200 MeV 
the (j) mean free path in hadronic matter is considerably 
smaller than what has been estimated so far. Finally, we 
study the implications of these findings for the cj) yields. 



II. THE HIDDEN LOCAL SYMMETRY 
LAGRANGIAN 

The HLS model provides a natural framework for de- 
scribing the interactions between pseudoscalar and vec- 
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tor mesons. It is based on the fact that a [U(3)lx 
U(3)r]/U(3)v non-Unear sigma model is gauge equivalent 
to another one with [U(3)lx U(3)R]giobai x [U(3)v]iocai 
symmetry. The most general Lagrangian possessing this 
symmetry and made with the smallest number of deriva- 
tives is |l6| 
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(2) 
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where /,r stands for the pion decay constant and a is an 
arbitrary parameter; it is usually set to a = 2, which al- 
lows to recover the standard vector meson dominance ex- 
pression (VMD) (l6[ |l^, although a slightly bigger value 
has been extracted from experiment p^ . The symbol 
denotes flavor trace. The covariant derivatives of the 
fields £_L.,R are 
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with being the dynamical gauge bosons of the hid- 
den local symmetry, further identified with the nonet of 
vector mesons 
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Working in the unitary gauge i.e. choosing S,l,r such 
that 



(6) 



£a becomes identical to the lowest order chiral La- 
grangian 



(7) 



with U — S,'^; $ is the octet of pseudoscalar Goldstone 
bosons 
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The physics contained in Cy becomes clear when ^ is 
expanded up to the term linear in $ 



(9) 



Then, 



Cy = iasV^^;^"" - d,^] , V^}) , (10) 

where [,] and {,} stand for commutator and anticom- 
mutator respectively. Vector mesons have acquired mass 
™y — 0.9^ fi by spontaneous breakdown of the hidden 
local symmetry. One now assumes that a kinetic term 
for them is generated by the underlying QCD dy namics 
or quantum effects at the composite level |7|. Thus, 
we have 



(11) 



being the nonabelian field strength tensor 

F^, = - d,V^ - ig [V;, K] . (12) 

Substituting Eq. into (|ll|) we obtain, apart from the 
standard kinetic terms, the interactions of 3 and 4 vector 
mesons 

Cvvv ^i9{{d^V,-d,V^)[V^,V^]) (13) 



Cvvvv = -\{[V^..VA[V^,V^]) (14) 

As we shall see, these terms generate large contribu- 
tions to the total cross section of the (j) meson with other 
hadrons, specially vector mesons. 

In order to account for deviation from the flavor sym- 
metry one should add [U(3)lX U(3)Rjgiobai breaking 
terms which do not affect the hidden symmetry. There is 
no unique way to do this. Different implementations are 
studied in Ref. . Here we adopt the scheme proposed 
in Ref. Ml , in which La and Cy are modifled as follows 



Ca + ACa ^-\f'{ {d^^l ■ el - D^^n ■ Cl 



l + ^LeA^l + ^ReA^D), (15) 



Cy + ACy = -\fl{ {D^a ■ el + D,^R ■ 

'l+aeyel,+?flevel)) , (16) 



with iA(v) — diag(0, 0, c^(,^)); c^(v) are real parameters 
to be determined. After fixing the gauge (Eq. ^) and ex- 
panding e (Eq. 1^) , one observes that the kinetic terms in 
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Ca have to be renormalized. This is achieved by rescal- 
ing @ 



(17) 



Such a redefinition leads to the foUowing relations for the 
pseudoscalar meson decay constants 



1 + U ■ 



Now, different vector mesons have different masses 
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1 + 2c,.. 



(18) 



(19) 



while the second term in Eq. ([10|), which describes the 
coupling of vector mesons to a pair of pseudoscalar ones, 
looks like 

Lvpp = '^^gm, df,^] , V^} (1 + 2ey)) . (20) 

Eq. ( p^ ) contains also a vector-vector-pseudoscalar- 
pseudoscalar contact term 

Cvvpp = -ag^ {V^V (ey$^ + 2$ey$ + <i>hv)) ■ 

Obviously, the flavor symmetry is also broken in the 
pseudoscalar sector. The corresponding Lagrangian can 
be written as 



(22) 

where x can be expressed in terms of tt and K masses 

X = diag(m2 , ml, 2m\ ~ ml) (23) 

assuming the same mass for both u and d quarks. 

It is known that the local chiral symmetry is broken 
at the quantum level. The anomalous part of the La- 
grangian in terms of effective degrees of freedom is ob- 
tained assuming that the anomaly at composite level 
should coincide with that at constituent level. In the 
framework of HLS, vector mesons were incorporated to 
the anomalous Lagrangian by Fujikawa et al. ^l|. It 
contains the vector-vector-pseudoscalar interaction. In 
Ref. |2^, the flavor breaking effect was included by in- 
troducing a term {£,l^wz£}ji + S.R^wzi\), with e.^z = 
diag(0, 0, Cvvz), so that the total anomalous Lagrangian 
reads 

Cvvp + l^Cvvp 

= 25^^pe^'^^" (a^K (1 + 2e„,) d^V^<^) . (24) 
The coupling constant is fixed by the anomaly 
35^ 



5v 



(25) 



and c„z is directly obtained from the ratio between the 
experimental decay widths of K*'^ K^j and K*^ 



III. COLLISION RATES AND MEAN FREE 
PATH 

A. General formulae 

The propagation of (p mesons through the hot hadronic 
gas depends on how often they interact in their way out 
of the fireball. For a given binary reaction 1 -f 2 ^ 
3 + 4 , the number of collisions per unit time at a given 
temperature is described by the average collision rate 
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(2^)V4)(pi +P2 -P3 -P4)/l/2(l +/3)(1 + h) (26) 



where \M.\ stands for the amplitude squared, summed 
over final spins and averaged over initial ones; di are the 
spin degeneracy factors di = 2Si + 1 (we take care of the 
isospin degeneracy by considering each isospin channel 
independently), while fi represent the momentum distri- 
butions. The particle number density is given by 



(27) 



Approximating all momentum distributions by Boltz- 
mann 



-Ei/T 



-Ei/T 



1 



(28) 



and, identifying the particle labeled 1 with the (j) meson, 
one gets 

d2 



.(2) 
" coll 



(T) 



ATT^mlK2{m^/T) 



ds^spl„Ms)KiiV~s/T) (29) 

with s and Pcm being the center-of-mass (CM) energy 
squared and three- momentum respectively and (t(s), the 
total cross section; index (2) denotes the particle that 
collides with the (j). 

The most straightforward way to estimate the rele- 
vance of rescatterings in the medium is to compute the 
mean free path and compare it with the size of the sys- 
tem. For the mean free path we use 



A(r) = 



r* 



(30) 



V is the average velocity of the (j) in the medium 



..iLffELj^M.IEE^^ (31, 

ni J {2ny El m^^K2{m^/T) 

the expression on the right is obtained assuming a Boltz- 

mann distribution. F'^^j — J2a ^''coii arises from the con- 
tribution of all relevant hadronic reactions. In the present 
calculation we have a = tt, K, p, lo, K* , cf). 
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B. Binary Reactions 

Within the HLS model described in Section II one can 
compute ah possible binary reactions of 4> mesons with 
TT^ K, p,LU, K* ^ (f) at tree level. The interaction vertices 
are given by the VPP, VVRVVV WW and WPP La- 
grangians defined in Eqs. ( po| , |2^ , |l3|, ^ respectively. 
There are three classes of Feynman diagrams, which are 
depicted in Fig. |^ and represent s-channel (s), t-channcl 
(t) and contact (c) reaction mechanisms. The reac- 




(s) 



K, K* 



(t) 



No. 


Reaction 




Channels 


1.1 






K + K 


t(K') 


1.2 






K + K* 


t{K,K'), c 


1.3 






K* + K* 


t{K,K') 


2.1 




K - 


n + K 




2.2 




K - 


^ p + K 


sMK,K*) 


2.3 




K - 


-*uj + K 




2.4 




K - 


^<j, + K 


s,t{K,K*), c 


2.5 




K - 


-^TV + K* 


s{K,K*), tiK*) 


2.6 


6 + 


K - 


+ K* 


s,t{K,K') 


2.7 


(h -A- 


X - 




sMK,K') 


2.8 




K - 


^<j, + K* 


s,t{K,K') 


3.1 




n — ^ 


K + K 


t{K,K'') 


3.2 






K + K* 


t{K,K*) 


3.3 






K* + K* 


t{K,K*), c 


4.1 


rh + 




K + K 


t{K,K'') 


4.2 


rh + 




K + K* 


t{K,K'-), c 


4.3 






K* + K* 


t{K,K*), c 


5.1 




K* 


Tl + K 




5.2 


fA 4- 


K* 


-^p + K 


sMK,K*) 


5.3 


rh 4- 


K* 


-^LU + K 


s,t{K,K') 


5.4 


rA 4- 


K* 


+ K 


s,t{K,K') 


5.5 


rA 4- 


K* 


-^TV + K* 


s,t{K,K') 


5.6 


ni 4- 


K* 


-^p + K* 


s,t(if, A"*), c 


5.7 


+ 


K* 


-^LU + K' 


s,t(iC, A"*), c 


5.8 


+ 


K* 


^<j, + K* 


s,t{K,K*), c 


6.1 


+ 


0- 


K + K 


t{K,K'), c 


6.2 


+ 


0- 


K + K* 


t{K,K''), c 


6.3 




0- 


K* + K* 


t{K,K*), c 



TABLE I: List of binary reactions of rjj's with pseudoscalar (tt, 
K) and vector (p, lu, K* , <j)) mesons, s, t and c correspond to 
the diagrams in Fig. ^ while the particles in brackets denote 
the intermediate mesons. 



FIG. 1: Generic Feynman diagrams for the different reaction 
mechanisms; (s) stands for s-channel, (t) for t-channel and (c) 
for contact. OZI suppressed diagrams have not been consid- 
ered. 



tions considered in the present study are listed in Table |. 
We have neglected OZI suppressed reaction mechanisms. 
Therefore, only s and t diagrams with a strange particle 
{K or K*) in the intermediate state have been consid- 
ered. This is also the reason why s-channel mechanisms 
are only present in the case of (j) collisions with strange 
mesons. Notice that there are two kinds of contact terms: 
the WPP vertex, which appears as a result of symmetry 
breaking and the WW vertex obtained from the kinetic 
part of the Lagrangian. 

Explicit expressions of the Lagrangian for the different 
vertices are written in Appendix They are derived 



by substituting the matrices of Eqs . in the generic 
Lagrangians of Eqs. (^, 13, 14, 2^, taking the trace 
and rescaling the kaon field as shown in Eq. (17^. From 
those Lagrangians it is straightforward to compute the 
Feynman rules for the vertices. We use the standard 
expressions for the meson propagators 



DK{q) 



q2 _ jj^2 



K 



(.1) = - g""-- 



'AT* 



(32) 



m 



K' 



tmK^Tp:\q^) 



(33) 



^tot) 



i stand for the total decay widths of the interme- 

diate K, K* and include all decay channels that are open 
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for a given value. Using the Lagrangians introduced 
above one gets 



a—p,uj,< 



+ rK^aK'iq'')0{q^ -ima+m},f)} (34) 



and 



^(tot) 



a — 7T,p.UJ,' 



^aK*iq^)0{q^~{ma + m*j,f)} , (35) 



where 9{x) is the standard step function. Exphcit ex- 
pressions for the different decay widths are hsted in Ap- 
pendix ^ By introducing such imaginary parts in the 
propagators, we are implementing an approximate uni- 
tarization for the s-channel diagrams. 

Since hadrons are extended objects, it is necessary to 
insert form factors that suppress high momentum trans- 
fers. This affects all t-channel Feynman diagrams. We 
adopt the widely used monopole form 



F.{q' 



A2 



A2 



K* 



q^ 



(36) 



assuming the same cutoff parameter A = 1.8 GeV for all 
species ^ . These form factors cause a strong reduction 
in the t-channel contributions to the collision rates. 

The amplitude squared for the WW term is quite 
large and increases very fast with s since it involves 



high powers of both g and s {g and 



As a conse- 



quence, higher order corrections become relevant as one 
goes away from threshold. We take approximately into 
account the resummation of s-channel loops in the con- 
tact amplitude for (j) + a—fb-\-c hy means of the sub- 
stitution 



ag 



[I + ag'^G{s, 7714,, ^a)]^^^ [1 + ag'^G{s, rub, rric)]^^^ 

(37) 

with 



G{s, rua, mb) 



.J_Pc 



8m^r 



(38) 

where Pcm is the momentum of the vector mesons a and b 
in the CM frame. The factor q = -1/2 for the (f>(j)K*K* 
vertex and a = 1/(2^2) for (f)pK*K* and (j)ujK*K*. 
This ansatz is inspired by the solution of the Bethe- 
Salpeter (BS) equation in the K-matrix approximation 
using only the part of the tree-level amplitude, which 
leads to an algebraic equation (see Appendix ^ for de- 
tails). G(s, m^, TO{,) vanishes at threshold, so that one 
recovers the tree-level amplitude, and becomes progres- 
sively larger with the increase of s, causing an effective 



suppression of the coupling constant g. At large s, our 
simple expression shall become inaccurate compared with 
the solution of the full BS equation with couple-channels 
effects incorporated, but one should bare in mind that 
high values of s are exponentially suppressed in the inte- 
gral that defines the collision rates in Eq. ( p9| ) . Therefore, 
this inaccuracy has a small numerical impact in the deter- 
mination of the collision rates, but is important to avoid 
artificially large contributions from cross section values 
above the unitarity limit. 

With these ingredients, one can compute the ampli- 
tudes; it is important to add them coherently since the 
interferences modify appreciably the total cross section. 
A caveat is in order regarding reactions 1.2 and 2.4. In 
both cases the exchanged kaon (in the t-channel) can 
be put on the mass shell, making the amplitude become 
singular. The singularity disappears if one takes into ac- 
count that kaons will develop an in- medium width 
However, studying how the modifications of propagators, 
vertices or masses and widths of the particles involved in 
the reactions influence the collision rates goes beyond the 
scopes of the present work. Hence, we just do not take 
them into consideration. In any case, their contribution 
will most likely shorten the mean free path even more. 



C. Results 

The parameters of the HLS model can be fixed using 
the experimental values of masses and pseudoscalar decay 
constants. We choose a=2 after VMD; then, using the 
p mass in Eq. ([l9|) and /^r = 92.4 MeV [|3| one gets 
g = 5.89. Once a and g are fixed, the remaining part of 
Eq. (|l^ ) gives Cv It ranges from 0.34 to 0.38 depending 
on weather determined from rriK' or m^. We adopt Cy = 
0.36. Next, from E q. ([18|) and using the experimental 
value fx/f-K = 1.22 |j7we obtain = 0.49. Finally, 
we take c^z = —0.1 as determined in Ref. |Q. Isospin 
degeneracy is assumed, i.e. the mass differences between 
mesons of the same species are not taken into account. 

The contributions of the different mesons to the 4> colli- 
sion rate as a function of temperature are shown in Fig. ||. 
For temperatures between 150 and 200 MeV, the collision 
rate is dominated by the K* followed by K and p, while 
the contributions from vr, uj and (j) are smaller. 

In order to understand why F^.^; ^ is bigger than F^.^] 

in spite of their mass difference let us recall that f|."|j can 
be cast as 



(39) 



Hence, the averaged rate factor (avrei) must overcome a 
ratio of (see Eq. |2^) 



riK mj^K2imK/T) 



(40) 



at T = 200 MeV this ratio is equal to 0.77. This number 
is much bigger than what one would naively expect. In 
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FIG. 2: Collision rates of <j) with tt, K, p, uj, K* and as a 
function of temperature. 



fact, the small number coming from the ratio of Bessel 
functions is almost compensated by the ratio of squared 
masses times the spin degeneracy. The former being due 
to the higher density of states for a heavy particle. There- 
fore, a larger r|,^j ^ can be achieved by a moderately larger 
cross section. The comparison between Figs. || and ^ 
show that the reactions involving all four vector mesons 
make the true difference between r|,^j ^ and r|.^) , and the 
mechanisms involving WW and VW vertices given in 
Eqs. (|l^, |l^ are responsible for this. 

The sum of all partial rates of Fig. ^ is shown in Fig. ||. 
It is considerably bigger than what has been predicted 
in all previous works 0, ||. In this case, the major 
novel ingredient is the dynamics involving different vector 
mesons as shown above. The dashed line represents the 
contribution of all ^-number changing reactions; that is, 
all reactions in Table | except the (quasi)elastic ones 2.4, 
2.8, 5.4 and 5.8. The inelastic reactions account for more 
than 80 % of the total collision rate. 

We have also studied the sensitivity of the result to 
the modification of parameters in the HLS model. It was 
found to depend appreciably on the value of g, which is 
not surprising since the cross sections involve high powers 
of this coupling (up to 17*). As mentioned in Section II, 
the value of a extracted from the experiment is slightly 
bigger than 2. If we take a — 2.4 then a lower 

limit of 5 = 5.38 is obtained. We fix the upper limit by 
using instead of rUp, getting g = 5.98. The region of 
possible values of Tcoii is bounded by the dotted lines in 
Fig. |. 

In Fig. ^ we show the cj) mean free path as a function of 
temperature. It goes below 2.4 fm at T > 170 MeV. This 
is much smaller than the typical size of the hadronic sys- 



FIG. 3: Contribution from the various reactions of <j!> with 
K* to the collision rates. The labels correspond to the reac- 
tion numbers in Table ^ The largest rate is from (j) + K* 
p -\- K* . Notice that the contributions of reactions 5.1 and 
5.3 as well as 5.2 and 5.7 almost coincide and are hardly dis- 
tinguishable in the plots. 
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FIG. 4: Contribution to the collision rate of </> with K from 
the different reactions listed in Table 0. 



tem created in heavy ion collisions (10-15 fm). Therefore, 
contrary to the common believe pi, the (/> mesons that 
are created after hadronization, at temperatures presum- 
ably between 170 and 200 MeV, will not leave the fireball 
without interacting; they will most likely be absorbed 
and reemitted in the hadronic phase. The uncertainties 
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FIG. 5: Total collision rate of (j) meson as a function of tem- 
perature (solid line). The dotted lines constrain the region 
of possible values when the model parameters are changed 
as described in the text. The dashed line is the contribution 
from (^-number changing processes alone. 



discussed above do not modify this conchision. 




150 155 160 165 170 175 180 185 190 195 200 
T [MeV] 

FIG. 6: Temperature dependence of the (j) mean free path in 
hot hadronic matter. Line styles have the same meaning as 
in Fig. |. 



IV. TIME EVOLUTION 

We now proceed to study qualitatively the time evolu- 
tion of the (j) density in the expanding hadronic fireball. 
It is governed by the following rate equation 



(41) 



where = 7(1, f) is the four velocity, defined in terms 
of the Lorentz factor 7 and the fluid velocity u\ stands 
for the source term 



a,b,c 



a^b,c 



[o'bc^4,aVbc)'nbnc 



tl>Vbc)nbnc , 



(42) 



b,c 



b,c 



which takes into account (/)-number changing processes: 
binary collisions, decay, recombination. In this approach, 
the evolution of momentum distribution is not investi- 
gated. Kinetic equilibrium (and the applicability of the 
hydrodynamical description) is assumed, but not chemi- 
cal equilibrium ||2^, |2^. 

In principle, Eq. ( pl[) is just one in a system of cou- 
pled equations that describes the time evolution of the 
densities of all species. Such a treatment would require 
the knowledge of the reaction cross sections for all the 
mesons of the model. We shall rather perform the simpli- 
fying assumption that, at hadronization temperature, all 
particles (tt, p, w, K*, (p) are in chemical equilibrium, 
and all of them except the <j) remain in it until freezeout. 
In chemical equilibrium, detailed balance holds, i.e. 



{(^<Pa^bcV^a)n;pnl^ = {(Tbc->4,aVbc)nl''n'"' 



and 



{T^^bc)n'^ = {(Tbc->4,Vbc)nl''n'"^ 



(43) 



(44) 



Substituting these equations in Eq. (^) and imposing 
that particle number densities for all species except 0's 
take their equilibrium values: 



eq 

na,b,c — n-a,b,c ' 



one obtains 



where n!^{T) comes directly from Eqs. (0^, |2q) 



eq 



(T) 



\ T J 



(45) 
(46) 

(47) 



r(T) denotes the number of interactions per miit time 
and is split in two parts 



r{T)^T,ouiT)+TdeciT). 



(48) 



^coU is the collision rate calculated in the previous Sec- 
tion. Notice however that the right hand side of Eq. (|4l) 
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does not hold for collisions of two 4> mesons </) + ^ 6 + c 
(reactions 6.1-6.3). The reason being that (f)'s are not in 
chemical equilibrium. For our simple estimate we thus 
ignore those processes which are very small (see short- 
dashed line in Fig |^). Analogously, the contribution of 
(quasi)elastic reactions <j) + a —t (p + c (2.4, 2.8, 5.4, 5.8) 
to the source term ^ vanishes if particles a and c are in 
chemical equilibrium. As we have already stated, these 
reactions account for less than 20 % of the total collision 
rate in Fig ||. 

^dec{T) is the average free decay width. For the sake 
of consistency we only consider decays into kaon pairs, 
which account for 83 % of the total width. The general 
expression for Tdec is a t rivial modification of Eq. (26). 
Proceeding as in Section III A one gets 



^dec{T) — Fq 



(49) 



where Fq = 3.7 MeV. Obviously, at high temperatures, 
Tdec is negligible compared to TcoU, but the later drops 
fast when the system cools down. 

In central relativistic heavy ion collisions the distribu- 
tion of matter is approximately uniform in rapidity, at 
least in the central rapidity region, and the geometry of 
the collision is cylindrically symmetric |^. Therefore, it 
is convenient to adopt as coordinates {r, ry, r, Lp\ where r 
and if are the transverse radius and polar angle, while r 
and rj represent the longitudinal proper time and space- 
time rapidity respectively 



I, t + z 

11 = — Ul 

' 2 t- z 



(50) 



The assumptions of a radial transverse expansion and 
Lorentz invariance in the central region imply that = 
u"^ = Q. Next, assuming a uniform density distribution 
in the transverse plane and averaging over the radial co- 
ordinate (analogously to what was done in Refs. j2^, |2^ 
for the spherically symmetric case) we get 



1 



d 



{TR^T)n^(u^)} = -F(T) (n^ 



(51) 

which is to be solved with the initial condition n^{To) = 
n'^^(To), tq and To being the hadronization time and tem- 
perature; R(t) is the transverse radius of the system. No- 
tice that tttR^ is nothing but the volume of the expand- 
ing fireball at mid-rapidity. The averaged r component 
of the four velocity is given by 



RHr) 



R{r) 



dr ru^ (r) . 



(52) 



Following Ref. 1^^, we assume that the flow vector 
can be constructed from two independent boosts, one 
in the longitudinal and another in the radial direction. 
Then, it is easy to see that at mid-rapidity 



1 



7r 



(53) 



A reasonable ansatz for radial velocity is |^ 
In the case of a constant (3s 



is a constant too. Using the boundary condition 



(54) 



(55) 



(56) 



one gets 



Rir) ^ (3siT - To) + Ro , (57) 
1.2 yl^/-^ fm is the radius of the colliding 



where _Ro 
nuclei. 

Finally, in order to solve the rate equation we need 
the relation between time and temperature, which can 
be derived from entropy conservation 

(su") = . (58) 

Repeating the same procedure as for one finds 

1 d 



tR'^{t) dr 



{TR^{T)n^{u^)] =0. 



(59) 



This equation can be easily solved, obtaining the follow- 
ing implicit expression for t(T) (or vice versa) 

tR\t){u-'){t)s{T) = toRI{u-){to)s{To) . (60) 

If {u^) does not depend on r, as for the ansatz consid- 
ered above, Eq. (|60| ) reflects that the total entropy in the 
central region of the collision is conserved. 

The functional dependence of the entropy upon the 
temperature depends on the properties of the system un- 
der consideration. We describe it as an ideal gas of tt, 
p, LO and K* obeying Boltzmann statistics 



s{T) 



^2 



(61) 



with i = {tt, K, p, oj, K*}; g denotes the degeneracy (both 
spin and isospin). The underlying idea is that the equa- 
tion of state of a gas of free hadrons and resonances 
should mimic the one of a dense hadron gas . 

Fig. 1^ shows the r dependence of the ratio of the cj) yield 
at mid-rapidity N{t) = TTTR'^(T)n^{T) and the number 
of them right after hadronization A^o = tqtt RqTi'^ (Tq) . 
We have used tq = 1 fm. To = 190 MeV, f3s = 0.6, a = 1 
and A = 197 (Au). 

At the early stages of the expansion, the system is hot 
and the collision rates are high, so that 0's are pushed 
towards equilibrium making its number decrease fast. As 
the system cools down, the collision rates become small 
and the contribution of the collision rates to the ratio 
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FIG. 7: Time dependence of the ratio of the <j!> number to its 
value at hadronization. The dashed hne shows the contribu- 
tion from inelastic collisions while the long-short dashed one 
accounts for free decay (and recombination) The combined ef- 
fect of both collisions and free decay is described by the solid 
line. 



starts to saturate. The effect of the decay is negligible 
below T = 2 fm, not only because Tdec <C Tcoii but also 
because detailed balance works more efficiently than at 
lower temperatures. We have taken a freezeout temper- 
ature of Tf = 100 MeV, which corresponds to a life time 
of Tf = 8.3 fm. If one assumes that only those 4> mesons 
present at freezeout are detectable via kaon pairs, then 
the (j) yield is reduced in 20 % with respect to hadroniza- 
tion. 

This model exhibits a notable dependence on the value 
of (3s, as illustrated in Fig |8[ The bigger Ps, the faster 
the temperature drops and, hence, the closer N/Nq gets 
to one. The case f3s = corresponds to the Bjorken 
limit, where the transverse expansion is neglected. For 
our choice of the equation of state, this is a rather unre- 
alistic situation because the system lives so long that the 
transverse expansion can not be neglected. For more re- 
alistic values of /3s we get a ratio N/Nq, at Tf = 100 MeV, 
between 0.78 and 0.86. 

We have also considered different hadronization tem- 
peratures, as shown in Fig. ^. If the hadronization tem- 
perature is low, the collision rates are small from the very 
beginning, and the interactions become less relevant. 

The present simple description of the time evolution of 
the (j) number in an expanding hadronic gas is intended to 
illustrate the consequences of the short (j) mean free path 
that we have obtained. A realistic analysis of the spectra 
would require a more sophisticated treatment. In partic- 
ular, the fact that the (j) collision rate is dominated by 
0-number changing reactions implies that mainly those 
processes are responsible for maintaining the kinetic equi- 
librium via detailed balance. Therefore, if the system is 
out of chemical equilibrium, detailed balance no longer 
holds and the system will also depart from kinetic equi- 
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FIG. 8: Ratio of the cf> number to its value at hadronization for 
different values of the transverse flow velocity at the surface. 
The end point of the lines correspond to a temperature of 
100 MeV. 
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FIG. 9: Same as Fig. ^ but for different hadronization tem- 
peratures. 



librium. This means that both (j) temperature and flow 
velocity will differ from those of other species with high 
rates of elastic scatterings like pions. 



A. The role of chemical potentials 

Up to now we have not introduced chemical poten- 
tials (fugacity factors) for any of the particle species 
under consideration. However, it is well known that 
when the number changing processes are not effective, 
particle numbers are fixed and chemical potentials asso- 
ciated with these conserved quantities appear 
For pions and kaons values of /ijr = 60 — 80 MeV and 
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IJ,K = 100 — 130 MeV are reached at freezeout tempera- 
tures between 110 and 120 MeV [||. 

In order to study the effect of such conserved me- 
son numbers in the (f> yield we assume that fj,i — at 
hadronization for all species and grows to a value fj,{ at 
freezeout. A simple linear interpolation is used in be- 
tween: 



Tq-T 
To-Tf 



f 



(62) 



We take fil = 70 MeV, = 115 MeV with Tf = 
115 MeV to be consistent with the values quoted above. 
Since p ^ tttt is a fast process, one can say that tt 
and p mesons are in relative chemical equilibrium i.e. 
p.p = 2p,T^ In the case of the w, a recent calcula- 

tion has obtained a large collision rate for the reaction 
LOTT TTTT at high temperatures | |3^ ; therefore, we take 
fJ-Lj = fJ-TT- With the K* the situation is uncertain because 
its decay width into Ktt (51 MeV) might or might not be 
large enough to ensure a relative chemical equilibrium. 
Since the K*^s are important for the (j) yield due to their 
large contribution to the collision rate, we consider here 
two possibilities. The first is the assumption of chemical 
equilibrium via K* ^ Kn, so that p.K-' = P-k + Ptt and 
the second is that the K* has its own chemical potential, 
which we keep equal to the initial one pK* = for the 
sake of simplicity. 

The presence of chemical potentials implies that the 
condition (EHI) should be replaced by 
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FIG. 10; Time dependence of the ratio N/Nq at zero (solid 
line) and nonzero (dashed line) chemical potentials, as de- 
scribed in the text. The dotted line stands for the case of 
PK* = while keeping 



therefore, a faster cooling and earlier freezeout. If one 
sets pk' to zero while keeping the previous values for the 
other mesons (dotted line) then the tendency is inverted, 
and direct reactions are the most efficient. However, the 
final (j) number is still bigger than in the pi = case due 
to the faster cooling of the hadron gas. 



fi en 



na.b,c — '"nhr — '•'a,b,c^ 

leading to the following rate equation 



(63) 



-Tdec [ n. 



l^a(T) „„ Mi,(T) + Hc(T) 



(64) 



The equation of state should also be modified. Now, 
instead of Eq. (^ij) , we have 



Pi{T) T dpi{T) 
m," m, dT 



K2[^)\. (65) 



The result for the ratio N/Nq is given in Fig. |To|. The 
solid line reproduces the calculation without chemical po- 
tentials. It corresponds to the solid line in Fig. ^but ends 
earlier due to the higher freezeout temperature consid- 
ered here. The introduction of chemical potentials causes 
an increase in the number of (/)'s as shown by the dashed 
line. In this case, the inverse reactions are favored keep- 
ing the (j) number higher. The higher density at a given 
temperature also translates into a higher pressure and. 



V. SUMMARY 

We have studied (f> interactions with a hot hadronic 
medium composed of tt, K, p, uj, K* and (p using the 
Hidden Local Symmetry model. In this way, we could 
take into account many vertices that are allowed by the 
symmetries of strong interactions but whose couplings 
cannot be directly determined experimentally. This is 
the case of three and four vector meson vertices which 
are responsible of the large collision rates of 0's with 
other vector mesons, specially K* . As a consequence, we 
have shown that the (p mean free path at temperatures 
above 170 MeV is between 2.4 and 1 fm i.e. much smaller 
than the typical size of the hadronic system created in 
relativistic heavy ion collisions. 

The implications of this result for the cf) yield has been 
investigated by solving the rate equation for the (j) den- 
sity assuming kinetic but not chemical equilibrium. The 
high collision rates at the early stages of the hadronic 
evolution tend to maintain the equilibrium, causing a re- 
duction of the (j) number with respect to hadronization. 
This decrease ranges from 5 % to 45 % depending on 
hadronization temperature, freezeout temperature and 
flow velocity. Finally, we have studied how the departure 
from chemical equilibrium of tt, p, K and K* mesons, 
taken into account the by introducing chemical poten- 
tials, influences the time evolution of the d> number. 
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APPENDIX A: FULL SET OF LAGRANGIANS 

The explicit expressions of the Lagrangians describing 
the interaction vertices used in this work are listed below. 



pKK 



* 4 1 + Ca 



(Al) 



C^KK = i^^— uj. (K-df'K+ - K+df'K- + K^df^K^ - K^d^K^) 
4 1 -I-C4 ^ ^ 



(A2) 



£^KK = i \+'2cv ^K+dt'K- - K-d^K+ + K'^dt'K^ - K'^d^K'') 



(A3) 



ag 1 -f Cv 



4 yrr^ 



From Eq. (|2|): 



+d^pl {K-dxK*+ + K+dxK*- - K^dxKf - K^dxKf) } , 



(A4) 



(A5) 



K'K 



9vvp 



2 vr+^ 



e^-'^^a.c^, (K-dxK*+ + K+dxK*- + K'^OxK*" + K°dxK*°) 



(A6) 



= gvj.l + 2cw.^^.,^ {K-dxK;+ + K+dxK- + K^dxK^ + K^dxK^) , 



(A7) 



9vvp 



(1 + 2c^^)e^'^^" {V2n+d^Kt'>dxK*- + V2tt- d,J<*+dxK, 



*o 



^n'^d^Kt+dxK*- - 7r'>d,,K;'>dxKf) . 



(A8) 
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From Eq. 

[{d^K*- - d,K*-)K*+- - {d^K*+ - d,K;+)K*-] - d^^p°''{K;-K*+ - K*+Kt-) 



C^K'K' = »| {^^ [{^^,Kt- ~ d,K*-)K*+- - {d^Kt+ - d,K*+)K*-] - d^^tu- {K;- K;+ - K*+Kl-) 



(K*^K*+ ~ K*+K*' 



From Eq. (|l|): 



K 



(A9) 



(AlO) 



(All) 



(A12) 



(A13) 



From Eq. 



C 



4>4>K*K* 



r 

' 2 



ag 



4v2 \/1 + Ca ^ 



(A14) 



(A15) 



(A16) 



APPENDIX B: K AND K* DECAY WIDTHS 

For the kaon: 
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APPENDIX C: BETHE-SALPETER EQUATION 
FOR THE WW AMPLITUDE 



^K-^pK'is) 



^K-.ujK'{s) 



167r 1 + c, 



PcmiSii^K'^rnp) , (B4) 



o2 1 

i/vvp ^ 3 



16tt I + Ca 



pU.s,mK',m^), (B5) 



i K^^K' (S) — 



in 1 + c. 



Pcm{s,mK',m4,) . 



(B6) 



For the K* 



a^g^ (1 + Cy)^ plmjs, mK,m^) 
^K'^7tk[s) = — — , (B7) 



TK^^.Kis) = f;^" pl„,is,mK,mp), (B8) 

iOTT i + Ca 



9vvP 1 

487r 1 + 



plmis,mK,m^) , (B9) 



9vvp (1 + 2Cvt,z)^ 3 



247r 1 + Ca 



Pcm.{s,mK,m^) , 



(BIO) 



Tif.^^if. (s) = -p7^(l + 2c„zfpl^^{s,mK',m^) , 



IGtt 



g^ pUs,mK.,m,) + iml, 



„3 



(Bll) 



Pcm{s,mK',mp) , (B12) 



5^ Pcm(s>'7lK*,TO„) + 3m^. 



247r 



Pcm ( 

s,mK',mu,) ,(B13) 



)i<--(s) 



127r 



plmis,mK',m^) . (B14) 



Here, s is the invariant mass squared of the K or K*, 
which coincides with the CM energy squared for s- 
channel diagrams, and pcm{s,mi,m2) is the CM three- 
momentum of the outgoing particles. 



Let us consider the elastic reaction (j){pi) + K*{p2) — > 
(f'iPs) + ^^*iPi)- The contact term amplitude, given by 

Mfuii = ag'^{ep{pi)e^{p2)ef,{p3)e^{p4) 

+ ef,{pi)ep{p2)e^{p3)£^{p4) 

- 2e^(pi)ej.(p2)e^(P3)eM(P4)} , (CI) 

where e^(p) is the polarization vector and a = —1/2, 
is too involved so that solving the BS equation for it 
becomes impracticable. Instead, we take a part of the 
full amplitude 

Mpart = ag'^ep{pi)ep(p2)e,j{p:i)ey{pi) (C2) 

for which the equation becomes algebraic. It is easy to 
see that the solution looks like 



'■part 



fis)£p{pi)ep{p2)e^{p3)e^{p4) 



with 
and 



ag- 



1 + ag^G{s, m,j„ rriK' ) 



(C3) 



(C4) 



(27r)4^^* 

Pi +P2 



Pi +P2 



(C5) 



Neglecting the real part of the propagators (K-matrix 
approximation) 



1 



(C6) 



one obtains Eq. ( pq ) for G. Then, we extrapolate Eq. (C4) 
to the case of inelastic reactions assuming that it can 
be separated into to factors: one corresponding to the 
incoming pair of particles and another to the outgoing 
one, arriving this way at Eq. (|3^). 

Even if AApart is not the correct tree-level amplitude, 
we believe that using Eq. (C4), or (|3^), instead of the 
factor ag^ takes fairly well into account the corrections to 
M. full arising from the resummation of s-channel loops. 
We have actually studied the non-relativistic case, where 
the BS equation for M fuii also becomes algebraic, finding 
that the ratios M'^^n/M^f^H and M^^^JM^^^t are very 
similar, at least for moderate values of s. 



14 



[1] J. Rafelski and B. MuUor, Phys. Rev. Lett. 48, 1066 

(1982). [19 
[2] A. Shor, Phys. Rev. Lett. 54, 1122 (1985). 
[3] M. Asakawa and C. M. Ko, Nucl. Phys. A572, 732 [20 

(1994). 

[4] C. Song, Phys. Lett. B388, 141 (1996), hep-ph/9603259. [21 
[5] M. Asakawa and C. M. Ko, Phys. Lett. B322, 33 (1994), 

nucl-th/9401016. [22 
[6] C. M. Ko and D. Seibert, Phys. Rev. C49, 2198 (1994), 

nucl-th/9312010. [23 
[7] K. HagUn, NucL Phys. A584, 719 (1995), nucl- 

th/9410028. [24 
[8] W. Smith and K. L. Haghn, Phys. Rev. C57, 1449 (1998), [25 

nucl-th/9710026. 

[9] S. V. Afanasiev et aL (NA49), Phys. Lett. B491, 59 [26 
(2000). [27 
0] M. C. Abreu et aL (NA38), Phys. Lett. B368, 239 (1996). [28^ 
1] M. C. Abreu et aL (NA50), J. Phys. G27, 405 (2001). 
2] D. Rohrich, ,L Phys. G27, 355 (2001). [29 
3] K. Safarik, ,L Phys. G27, 579 (2001). 
4] S. C. Johnson, B. V. Jacak, and A. Drees, Eur. Phys. J. [30 

C18, 645 (2001), nucl-th/9909075. 
5] M. Bando, T. Kugo, and K. Yamawaki, Phys. Rept. 164, [31 
217 (1988). 

6] M. Bando, T. Kugo, and K. Yamawaki, Nucl. Phys. [32 

B259, 493 (1985). 
7] M. Bando, T. Kugo, S. Uehara, K. Yamawaki, and [33 

T. Yanagida, Phys. Rev. Lett. 54, 1215 (1985). 
B] M. Benayoun et aL, Eur. Phys. J. C2, 269 (1998), hep- 



ph/9707509. 

M. Benayoun and H. B. O'Connell, Phys. Rev. D58, 
074006 (1998), hep-ph/9804391. 

A. Bramon, A. Grau, and G. Pancheri, Phys. Lett. B345, 

263 (1995), hcp-ph/9411269. 

T. Fujiwara, T. Kugo, H. Tcrao, S. Uehara, and K. Ya- 
mawaki, Prog. Theor. Phys. 73, 926 (1985). 
A. Bramon, A. Grau, and G. Pancheri, Phys. Lett. B344, 
240 (1995). 

C. Caso et aL (Particle Data Group), Eur. Phys. J. C3, 
1 (1998). 

H. Leutwyler and M. Roos, Z. Phys. C25, 91 (1984). 
T. Biro, H. W. Barz, B. Lukacs, and J. Zimanyi, Phys. 
Rev. C27, 2695 (1983). 

C. M. Ko and L.-H. Xia, Phys. Rev. C38, 179 (1988). 

J. D. Bjorken, Phys. Rev. D27, 140 (1983). 

E. Schnedermann and U. Heinz, Phys. Rev. C47, 1738 

(1993). 

J. Sollfrank et aL, Phys. Rev. C55, 392 (1997), nucl- 

th/9607029. 

H. Bobic, P. Gcrbcr, J. L. Goity, and H. Leutwyler, NucL 
Phys. B378, 95 (1992). 

C. Song and V. Koch, Phys. Rev. C55, 3026 (1997), nucL 
th/9611034. 

R. Rapp and E. V. Shuryak, Phys. Rev. Lett. 86, 2980 
(2001), hep-ph/0008326. 

R. A. Schneider and W. Weise, Phys. Lett. B515, 89 
(2001), hep-ph/0102189. 



